Gluon Condensate in Pion Superfluid beyond Mean Field Approximation 
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We study gluon condensate in a pion superfluid, through calculating the equation of state of 
the system in the Nambu-Jona-Lasinio model. While in mean field approximation the growing pion 
condensate leads to an increasing gluon condensate, meson fluctuations reduce the gluon condensate 
and the broken scalar symmetry can be smoothly restored at finite isospin density. 
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Quarks and gluons condense in the vacuum of Quan- 
tum Chromodynamics (QCD). From lattice QCD calcu- 
lations and effective QCD models in hot medium, it is 
widely accepted that the quark condensate (V'V') which 
is the order parameter of the chiral symmetry restoration 
drops down at finite temperature. The gluon condensate 
(GJ^^G'{^'^)P, @| which describes the degree of the scale 
symmetry breaking is, however, not so optimistic. 

The gluon condensate at finite temperature is inves- 
tigated in instanton modelQ, renormalization groupQ, 
QCD sum rule@, and effective QCD models at low 
energy [^-[l^. While the results in these calculations are 
quantitatively different, they show the same temperature 
trend of the gluon condensate: It is almost invariable at 
low temperature and starts to decrease rapidly around 
the critical temperature of QCD phase transitions. 

At finite isospin density, both the Lee-Huang- Yang 
model [H for a dilute Boson gas and the Nambu-Jona- 
Lasinio (NJL) modcl[14] show a surprising mean field 
result [l^: In the pion superfluid the gluon condensate 
drops down slightly only at very low isospin density but 
goes up and even exceeds its vacuum value when the den- 
sity is high enough. This result is qualitatively in agree- 
ment with the calculations for 2-color baryon matter and 
3-color isospin matter [T7| . A natural question is if 
this conclusion is still true when we go beyond the mean 
field. In this paper, we study the gluon condensate in 
pion superfluid in the NJL model beyond mean field. 

Neglecting the current quark mass m, the QCD La- 
grangian is invariable under the scale transformation 
il}{x) — >■ A^/^'0(Aa:) for the quark field and A^^{x) — >■ 
\A^{\x) for the gauge field. At classical level, the trace 
of the corresponding Noether current is d^J^^ — — 

mtptp. At quantum level, the running coupling constant 
as leads to a so-called anomaly term, the trace of the 
ensemble average of Tlf becomes exactly the trace of the 
energy-momentum tensor of the system, and therefore 
the matter parts of the quark and gluon condensates at 
finite temperature T and chemical potential are related 
to the energy density e and pressure p of the svstem[l8l|. 
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This relation tells us that the QCD condensates are con- 
trolled by the bulk properties of the system. Since it 
is difficult to directly calculate the QCD thermodynam- 
ics in non-perturbative region, this relation gives a way 



to qualitatively estimate the gluon condensate in effec- 
tive models at low energy where partons are not explicit 
constituents, if the model can reasonably describe the 
QCD thermodynamics. For instance, the gluon conden- 
sate has been investigated in nuclear matter (7Hlll] and 
in isospin matter [isl Hal with low-energy models. When 
we neglect the current quark mass m, the gluon conden- 
sate decouples from the quark condensate and is purely 
controlled by the thermodynamics of the system. While 
the gluon condensate for an ideal gas with e — 3p = is 
medium independent, it will be significantly changed for 
a strongly coupled system. From the lattice simulation 
at finite temperature |19|], the QCD system is a strongly 
coupled matter around the phase transition temperature 
Tc with e — Bp ^ 0. This is the reason why the gluon 
condensate drops down dramatically around Tc- 

The NJL model at quark leveljl^ has been successfully 
used to study chiral symmetry restoration, color super- 
conductivity and pion superfiuidity at moderate temper- 
ature and density. The flavor SU(2) NJL model is defined 
through the Lagrangian density 

r 2 2" 

£ = V' (^7^9^ - m + ^7o) + G (i/;?/;) + (V'^Ts'^V') 

"(2) 

where the quark chemical potential matrix /i = 

diag{^u, Md) = diag{^B/S + — P^i/'^) and the 

Pauli matrices r = (ti, T2, T3) are defined in flavor space, 
/is and /i/ are baryon and isospin chemical potentials, 
and G is the four-fermion coupling constant. The NJL 
thermodynamic potential can be separated into a mean 
field part and a fiuctuation part. 
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The mean field part flMF contains the mean field po- 
tential and the contribution from the a uasi-g uarks [2 Ij . 
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where the chiral condensate a = (V'V') and pion con- 
densate TT = V2{il'ij5T+4') with r+ = (ti -|- iT2)/y/2 are 



2 



determined by minimizing the potential, 
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and — E^ =p fJ'B/'i are the quasi-quark energies with 
E^ = 



{Ek ±^iI/2y + 4G^n^, Ek = Jk^ + and 



da 



dynamical quark mass Mq — m — 2Ga. 

In the NJL model, the meson modes are regarded as 
quantum fluctuations above the mean field. The two 
quark scattering via meson exchange can be effectively 
expressed in terms of quark bubble summation in ran- 
dom phase approximation [2^. In normal phase without 
pion condensation, the bubble summation selects its spe- 
cific isospin channel by choosing at each stage the same 
proper polarization, and the meson masses Mm (jn = 
(7, 7r_|_ , 7r_ , ttq ) which are determined by poles of the me- 
son propagators, 1 — 2GIlm7n{Mm,0) — 0, are related 
only to their own polarization functions Ilmm{<lo,(l)- In 
pole approximation, the meson contribution to the ther- 
modynamic potential can be expressed asj22| 
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with meson energies Em — -^/q^ -I- and meson isospin 
chemical potential /i7r± — ^^J'I and /z^g = /io- = 0. 

In the pion superfiuid phase, the quark propagator 
contains off-diagonal elements in flavor space, we must 
consider all possible isospin channels in the bubble sum- 
mation. In this case, all the possible polarizations form 
a matrix 11 in the four-dimensional meson isospin space 
with off-diagonal elements Hmn- While there is no mixing 
between ttq and other mesons, 117^0- = n7ro7r+ 



0, the other three mesons are coupled to each other. The 
explicit T, /i^ and /i/ dependence of all polarization ele- 
ments Hmn can be found in Appendix B of [2l|. When 
the system goes through the phase transition line and en- 
ters the normal phase, all the off-diagonal elements dis- 
appear automatically. 

The masses of the eigen modes of the Hamiltonian H 
in the pion superfiuid are defined through the poles of the 
meson propagator, det(l — 2Gll(Mg,0)) = which can 
be separated into 1 — 2Gn7ro7ro (M^^, 0) = for = ttq 
and det(l — 2GIl{Mg,0)) = in the three-dimensional 
isospin subspace for 9 = a, 7f_|_ , tt _ . Different from the 
normal phase where the meson modes cr, 7r+ , 7r_ , tto are 
eigen states of both the Hamiltonian T-L and the isospin 
operator 1^ — 1/2 J (Px^p^oTstp of the system, only Wq 

is still the eigen state of I3 (we still label it ttq in the 
following) , but ct, 7f + , 7f_ have no longer definite isospin 
quantum number. The eigen states of I3 are only re- 
lated to the diagonal elements n„i,„ and their masses are 
defined by 1 - 2Gn,„™(M„, 0) = 0. 

After taking bubble summation and Matsubara fre- 
quency summation, the fluctuation part of the thermo- 



dynamic potential can be generally written as [22 
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where the two polarization matrices are respectively de- 
fined in the top and bottom complex meson energy plane. 
An often used simplification to calculate £7^^ is the pole 
approximation, namely neglecting the scattering phase 
shifts and considering only the contribution from the 
quasi particles, like ^ for the normal phase. In this case, 
we have ri^L = X^e^e- To explicitly show the isospin 
dependence, we further make a transformation [2^ from 
the basis (ct, 7f-|_ , 7f _ ) to the basis (cr, -k^ , tt- ) . The el- 
emental states in the former basis do not carry definite 
isospin quantum numbers, but the later is constructed by 
the eigen states of the isospin operator 73. Since the two 
spaces are both complete, such a transformation will not 
lose any information. Taking into account the orthogonal 
condition for the two spaces, Q.0 can be expanded as a 
linear combination of VLm- Finally, we have 
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with the coefficients 

Cm^Y.\^e\, 
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where is a matrix defined in the three dimensional 
meson isospin subspace, 



— det(l - 2Gn(Afe,0)) 
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It is easy to see the normalization condition for the 
coefficients, J^m'^'m ~ J2e ~ '^^ means that only two 
of the three coefficients are independent. The coefficients 
Cm as functions of temperature at fixed chemical poten- 
tials are shown in FigHJ Their strong deviation from unit 
indicate a strong mixing of tr, 7r+ , 7r_ in the pion super- 
fiuid. For the Goldstone mode 7f_|_, its linear combination 
is |7f+) ~ l/A/2(|7r+) — |7r_))[23], and the two fractions are 
equal and medium independent. Therefore, at the criti- 
cal point the coefficient Cj^^ jumps up from 0.5 to 1 and 
C7r_ drops down from 1.5 to 1. For T > Tc in the normal 
phase, all the three coefficients are unit. For /i/ — 200 
MeV and = 600 MeV in FiglH Tc is about 110 MeV. 
It is necessary to note that the discontinuity of the coeffi- 
cients C7r_|_ and happens on the whole phase transition 
border. However, when we approach to the border from 
the pion superfiuid side, the pion condensate goes to zero 
continuously, and this can smooth the thermodynamics 
on the border, see the calculations below. 

Now we use the trace anomaly relation ([T]) to calcu- 
late the giuon condensate, under the assumption that the 
NJL model can describe reasonably well the QCD ther- 
modynamics in the pion superfiuid. From the thermody- 
namic potential relative to the vacuum Q.{T, fiB: fJ-i) — 
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FIG. 1: The linear coefRcients Cm for the transformation from 
the eigen states of Hamiltonian to the eigen states of isospin. 



ri(r, /XB, (ii) — fi(0, 0, 0), we obtain the pressure p = — 
and energy density e = —p + Ts + ^3^3 + IJ^i^i with the 
entropy density s = —dQ/dT, baryon number density 
n3 = —dfl/dfj,3 and isospin number density —dfl/d^j. 

Before we make numerical calculations, we first de- 
termine the parameters in the model. Since the NJL 
model is non-renormalizable, we can employ a hard three 
momentum cutoff A to regularize the gap equations for 
quarks and pole equations for mesons. In the following 
numerical calculations, we take the current quark mass 
7710 = 5 MeV, the coupling constant G = 4.93 GeV^^ 
and the cutoff A = 653 MeV. This group of parameters 
corresponds to the pion mass = 134 MeV, the pion 
decay constant f-^ — 93 MeV and the effective quark 
mass Mq = 310 MeV in vacuum. 

We show in Figl2]thc ratios for gluon, chiral and pion 
condensates, Rg = (s^G^^Gr)/(^G^.Gr)o, i?. = 
ct/cto and i?^ = tt/ctoj where (^G°yG^'')o and ctq are 
the condensates in vacuum, and {^G^j^^G^") , a and n are 
the total condensates including the vacuum and matter 
parts. To reduce the model dependence and focus on the 
medium effect, we take an empirical value for the vac- 
uum part of the gluon condensate, (^G°^G(^'')o=(360 
MeV)'*!!^] (the value of the vacuum part will not change 
the trend of the gluon condensate in the medium). At 
T = fiB = in the top panel of Fig 121 the ratio 
Jl^iP+FL^ calculated with the total thermodynamic po- 
tential 51 = flMF + ^FL, is a constant in the normal 
phase with /i/ < to^ and drops down monotonously in 
the pion superfluid phase with fij > m^. Therefore, the 
behavior of the gluon condensate at finite isospin density 
is qualitatively the same as in the case at finite temper- 
ature: The broken scale symmetry of the system is grad- 



ually restored in hot and dense medium. However, in 
mean field approximation, the gluon condensate behaves 
very differently. The ratio decreases slightly only 

in the beginning of the pion superfluid and then goes up 
monotonously and even exceeds the vacuum value when 
fii is high enough. For T = 50 MeV and = 600 MeV 
shown in the bottom panel of FiglJl while the mean field 
calculation is changed slightly, the finite temperature and 
baryon chemical potential effect results in stronger meson 
fluctuations, and the ratio R^^^^^^ drops down much 
faster. 
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FIG. 2: The scaled gluon condensates Rg'^ in mean field 
approximation and Rg'^'^^^ including quantum fluctuations. 
As a comparison, we showed also the scaled chiral condensate 
Ra and pion condensate R-^ at mean field level. The top panel 
and bottom panel correspond respectively to T — fiB = and 
T = 50 MeV, ^is = 600 MeV. 

The mean field result can be understood by the compe- 
tition between the chiral and pion condensates. At mean 
field level, the NJL Lagrangian density can be written as 



C 



MF 



= ip (17'' 9^ — m + /i7o) -I- 2G [aipip + Tr^'HsnV'] 
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- G(a2+7r2), (11) 

and the corresponding trace of the Noether cur- 
rent for the scalar transformation is = mtp'ip — 
2G (^aipip + 7r'(/'i75Ti'0) . Taking the identification of the 
trace of the energy-momentum tensor in QCD and in the 
NJL model, the gluon condensate is characterized only 
by the two condensates, 

= ^-^Ga^ {Rl + Rl ~ I) . (12) 

In the pion superfluid phase, the two ratios R^ and i?^ 
behave in an opposite way, Ra- drops down but Rtt goes 
up, and the trend of the gluon condensate is controlled by 
the competition between the chiral and pion condensates. 
When fi] is above but close to the critical point /ij — m^r, 
the chiral and pion condensates are equally important 
and their competition may result in a possible decreasing 
gluon condensate. However, when ^/ is large enough, the 
chiral condensate becomes small and the pion condensate 
dominates the system. In this case, the gluon condensate 
increases with increasing pion condensate. 

It is necessary to emphasize again that the trace 
anomaly relation ([T]) between the gluon condensate and 
the thermodynamics of the system is valid only at quan- 
tum level. At classical or mean field level, the relation is 
not true, and the scale symmetry of QCD is only explic- 



itly broken by the current quark mass m, (T^) = m{'ijj'ip). 
In the NJL model, the quantum fluctuations or the meson 
modes can not be neglected. At mean field level, there 
are only quarks in the model which control the thermo- 
dynamics only at high temperature and density. At mod- 
erate temperature and density around the chiral and pion 
superfluid phase transitions, both quarks and mesons are 
important. At low temperature and density, mesons be- 
come the dominant contribution to the thermodynamics. 
Therefore, we need quantum fluctuations to describe the 
system in the whole temperature and density region. 

In summary, we have studied the gluon condensate be- 
yond mean field approximation in a pion superfluid de- 
scribed by the NJL model. Since the trace anomaly rela- 
tion is valid only at quantum level, the quantum fluctua- 
tions in the model must be considered in the calculation 
of gluon condensate. At classical or mean field level, the 
growing pion condensate in the superfluid leads to a sur- 
prising increase of the gluon condensate. However, when 
the quantum fluctuations are included, the meson contri- 
bution dominates the thermodynamics of the system at 
low and intermediate temperature and density, and the 
gluon condensate becomes to decrease gradually in the 
pion superfluid. Therefore, the scale symmetry can be 
restored at both finite temperature and density. 
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